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ABSTRACT 


This Final Progress Report for NCC-1-178 presents the details of the 
engineering developement of an analytical/computational solution to the 
heavy ion transport equation in terms of a multi-layer Green's function 
formalism as applied to the Small Spacecraft Technology Initiative (SSTI) 
program. 

The mathematical developments are recasted into a series of efficient 
computer codes for space applications. The efficency of applied algorithms 
is accomplished by a nonperturbative technique of extending the Green's 
function over the solution domain. The codes may also be applied to the 
accelerator boundary conditions to allow code validation in laboratory 
experiments . 

Correltions with experiments for the isotopic version of the code with 
59 and 80 isotopes present for a two layers target material in water 
has been verified. 



INTRODUCTION 


Future NASA missions will be limited by exposure to space radiations 
unless adequate shielding is provided to protect men and equipments 
from such radiations. Adequate methods required to estimate the damage 
caused by such radiations behind various shields can be evaluated prior 
to commitment to such missions. 

From the inception of the Langley Research Center heavy ion (HZE) 
shielding program (refs. 2-4), there has been a continued, close relation- 
ship between code development and laboratory experiment (ref. 4). Indeed, 
the current research goal is to provide computationally efficient high 
charge and energy ion (HZE) transport codes which can be validated with 
laboratory experiments and subsquently applied to space engineering desigr 
In practice, two streams of code development have prevailed due to the 
strong energy dependence of necessary atomic/molecular cross sections anc 
the near singular nature of the laboratory beam boundary conditions (refs. 
5-7). The atomic/molecular cross section dependence is adequately dealt 
with by using the methods of Wilson and Lamkin (ref. 8), allowing effici- 
ent numerical procedures to be developed for space radiations (refs. 
7,9-11). Although these codes could conceivably be applied to the labo- 
ratory validation, methods to control truncation and discretization errors 
would bear little resemblance to the space radiation codes attempting to t 
validated. Clearly, a radical reorientation is required to achieve the 
validation goals of the current NASA space radiation shielding program, 
and such an approach in terms of a multi-layered Green's function formalis 
as an extension to a previously developed single-layered algorithm is the 
main thrust of this continuing research and is briefly described below. 


GREEN* S FUNCTION METHODS 


Transport Equations and Conventional Approach 

The transport equation for high energy heavy ions is usually simplified 
by assuming the straightahead approximation and neglecting the target 
secondary fragments (ref. 3 et. al) and is written as 
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where^(x, E)is the ion flux at x with energy E (MeV/amu) f Sj(E) is the 
change in E per unit distance, oj the total macroscopic absorption cross 
section and aj k the macroscopic cross section for collision of ion type 
k to produce an ion of type j . The solution to equation ( 1 ) is to be 
found subject to the boundary condition 


^>( 0 , E) = fj{E) 


which for laboratory beams has only one value of j for which fj(E) is 
not zero and that /j(E) is described by a mean energy spread a such that 

,j(E) = eXP [ -(E - B ° )2/7 ff2 ] (3 

The solution to equation (1) is given by superposition of Greens 
function Gjfc as 

B) = 

where Green function is a solution of 
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subject to the boundary condition 
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The above equations can be simplified by transforming the energy into 
the residual range as 
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and defining new field functions as 
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equation ( 5 ) becomes 
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with the boundary condition 
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and with the solution to the ion fields given by 
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Note*') is the range scale factor asujfj = i/ m r m and is taken as vj = Zj/Aj 
The solution to equation (11) is written as a perturbation series 
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fori/ m > Vj. If Uj > t'mas can happen in neutron removal, the negative of 
equation (16) is used and upper and lower limits of equation (17) are 
switched. The higher order terms are approximated as 
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Note that Qj m (x, rj, r m ) is purely dependent on x for i > 0, which we represent as Q^^x). (See 
ref. 3.) In terms of the above, the solution to equation (1) becomes (from ref. 3) 

( x > r j) = e ~ ajX fj ( r i + x) 
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In equation (21), r mu and r' jni are given by the upper and lower limits of the inequality of 
equation (17). The symbol F'mi^rn) rsfers to the integral spectrum 
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We now introduce nonperturbative terms for the summation in equation (21). 

* First > we reca11 that the 9 function of n arguments was generated by the perturbation solution 
of the transport equation neglecting ionization energy loss (ref. 1) given by 
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for which the solution is 
0jm(±) = 6 jm g{m) +-0j m g{j, m) + • • • 

It is also true that 
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for any positive values of x and y. Equation (26) may be used to 
propagate the function gj m {x) over the solution space from very thin 
shield solutions. Equation (14) is then rewritten as 
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and the approximate solution of equation (1) is given by 
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Green's Function Methods in a Shielded Medium 
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The major simplification in the Green's function method results from 



the fact that the scaled spectral distribution of secondary ions to a 
first approximation depends only on the depth of target penetration 
as seen in equations (16), (18), and (27). The first approach to a 
multi-layered Green's function will rely on this observation and assume 
its validity for multi-layered shields. 

Consider a domain labeled as 1 which is shielded by a second domain 
labeled as 2, the number of type j ions at depth x in 1 due to type 
m ions incident on domain 2 of thickness y is 

9x V J**v) = £ $,,*(*) (29) 

The leading term in equation (29) is the penetrating primaries as 
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where all higher order terms are in the bracket of equation (30). 
The first term of scaled Green's function is then 
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where p is the range scale factor for the two media 


P = R\j{E)/R 2 j {B) 

A sigle value for p is taken corresponding to 600 MeV/amu. The 


secondary contribution is similarly found by noting the equation (17) 
becomes 
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from which the average spectrum is evaluated. The full approximate 
Green's function is then 
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Equation (34) is the first approximation to the Green's function in a 
shielded medium (two layers) and can be modified to multiple layers. 
The first spectral modification is considered next. It is easy to show 
that the first collision term has the properties 
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These properties are used to correct the average spectrum as 



(1) , . 
v ± g|2jm(^y) 

Wm ~ Vj\(x + py) 


"f* b jm( x > y)( r m r m ) 


.( 1 ) 
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the first collision term of equation (34) and 
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is the midpoint of f' m between its limits given by equation (33). The 
term of equation (36) has the property that 
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ensuring that the first term of equation (36) is indeed the average 
spectrum as required. The spectral slope parameter is found to be 
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A similarly simple spectral correction could be made to the higher 
order terms. The spectral correction given in equation (39) is included 
in the present Green's function code. 

Solution for Laboratory Beams 

Using the boundary condition appropriate for laboratory beams given 
by equation (3), cumulative spectrum given by 

F i (E) = i[l-er/(^)] <40 


and the cumulative energy moment needed to evaluate the spectral 
correction given by 
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The average energy on any subinterval (E\,E 2 ) then 
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and the beam generated flux becomes 
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where® is evaluated using equation (42) with Ej, and £; 2 as the lower 
and upper limit associated with r’ m ^ and ^mxt • 


DISCUSSION OF RESULTS 

A series of algorithms in support of Small Spacecraft Technology 
Inititive (SSTI) based on the multi-layered Green's function 
formalism were developed and tested during this research project. 

In order to implement the theoritical formulations presented in the 
theory section of this research grant for multi-layered targets, the 
following tasks were undertaken: 

1. A direct conversion from formulation to coding based on the des- 
cribed procedures was performed. This resulted in a first gene- 
ration multi-layered computer codes with acceptable computational 
efficiency for the low order terms. 

2. Within the scope of the grant. Further improvments were carried out 
to make the higher order terms, as computational lly efficient as 
possible. The current version of these codes are computationally 
efficient for nearly all terms. 

3. The resulting coded package were further tested against the results 
of a recent two-layered experiments for correlation purposes and 
code efficency studies. 

4. The tested version of the two-layered algorithm was used as a 

mean of evaluating the shielding properties of different two-layered 
target materials as they are exposed to galactic cosmic radiations 


(GCR) . 
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